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Abstract 

In order to describe the fusion of two very heavy nuclei at near barrier en- 
ergies, a generalized Langevin approach is proposed, which incorporates the 
quantum statistical fluctuations in accordance with the fluctuation and dissi- 
pation theorem. It is illustrated that the quantum statistical effects introduce 
an enhancement of the formation of compound nucleus, though the quantum 
enhancement is somewhat less pronounced as indicated in the previous calcu- 
lations. 
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I. INTRODUCTION 

In recent years, there has been a great deal of interest in the synthesis of superheavy 
elements by means of heavy-ion fusion at near barrier energies [1]. Due to very low pro- 
duction probabilities, investigation of heavy-ion fusion reactions remains a challenging task 
both experimentally and theoretically. Theoretical investigations are mostly based on trans- 
port theory in which heavy-ion fusion is viewed as a diffusion process [2-5]. During the 
approach phase, a part of the kinetic energy of the relative motion dissipates into intrinsic 
degrees of freedom, colliding ions overcome the Coulomb barrier and a sticking configuration 
is formed. Subsequently, the system evolves inward over the conditional saddle to form a 
spherical compound nucleus. Heavy-ion fusion is described in terms of a few relevant vari- 
ables, which evolve according to a Langevin dynamics as in a typical diffusion process. In 
most of these investigations quantum statistical effects are ignored and a classical treat- 
ment is employed, in which friction and diffusion properties are related through classical 
fluctuation-dissipation theorem. Since superheavy elements are stabilized by the shell cor- 
rection energy, they should be synthesized at reasonably low energies, which corresponds to 
nuclear temperatures of the order of T = 0.5 — 1.0 MeV. On the other hand, the curvature 
of the conditional saddle is also of the order of HQ = 1.0 MeV. As a result, it is expected 
that the quantum statistical effects play an important role in the formation of a compound 
nucleus by diffusion along the conditional potential barrier [6-8]. In a recent work, using 
the density matrix formalism, we derived a generalized Fokker-Planck equation for the dis- 
tribution function of relevant collective variables with non-Markovian transport coefficients 
and we illustrated the quantum statistical effects on the formation of compound nucleus at 
low temperatures [9,10]. Here, we consider the same problem, but follow a different descrip- 
tion based on generalized Langevin approach. By projection on the collective space, it is 
possible to deduce stochastic equations of motion for collective variables. The friction term 
involves memory effect and the random force acts as a correlated noise and they are related 
according to quantum fluctuation-dissipation theorem [11,12]. As a result, the quantum 
statistical fluctuations are incorporated into the description and the approach is valid at 
all temperatures. In principle, the Langevin approach presented here is equivalent to the 
density matrix formalism, however it has certain advantages as compared to the description 
by the Fokker-Planck equation [13]. One thing is that, in realistic situations, numerical 
simulations of Langevin equation require much less numerical effort than the solution of the 
Fokker-Planck equation. Furthermore, diffusion coefficients presented here are modified by 
the friction mechanism, which was not considered in the density matrix approach [9,10]. 
In section 2, we derive generalized Langevin equations for the relevant collective variables. 
In section 3, we present an analysis of the Langevin equation. In section 4, in order to 
illustrate the effect of quantum statistical fluctuations, we present calculations to describe 
the heavy-ion fusion reactions at low temperatures. Finally, we give conclusions in section 
5. 

II. GENERALIZED LANGEVIN EQUATIONS 

It is possible to derive transport equation for the entire dynamics of the heavy-ion fusion 
process, starting from the entrance channel until formation of a compound nucleus. Here, 
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we want to investigate the influence of quantum statistical fluctuations on the formation 
probability of compound nucleus at near barrier energies. For this purpose, we consider a 
part of the dynamics, namely, evolution of the system along the conditional saddle from 
sticking configuration until formation of compound nucleus. For simplicity, we describe the 
evolution in terms of a single collective variable, namely the elongation variable, q, which 
approximately corresponds to the relative distance between the colliding ions. We consider 
a model, in which evolution of the system is described by a Hamiltonian of the form, 

H — H + H cot i + V coup (1) 

where H , H coU and V coup represent the Hamiltonian of the intrinsic nucleonic degrees of free- 
dom, the Hamiltonian of the collective variable and the coupling interaction of the collective 
variable with the intrinsic nucleonic degrees of freedom, respectively. For simplicity, we 
consider a harmonic form for the collective Hamiltonian H co u = p 2 /2M ± MQ 2 q 2 /2, where 
M is the mass parameter of the collective variable and MQ 2 denotes the magnitude of the 
curvature parameter of the potential energy, positive sign stands for a parabolic potential 
well and negative sign for a parabolic potential barrier. Furthermore, we assume that the 
coupling Hamiltonian has a linear form, V coup = qF. The classical equation of motion for the 
collective variable can be deduced from dp/dt = — (dH/dq), where (..) denotes an average 
over the intrinsic degrees of freedom, to give, 

j t Q(t) = ± P (t) and |p(t) ± MQ 2 q(t) = -Tr(Fp) . (2) 

Here, the quantity on the right hand side denotes the force of the intrinsic degrees of freedom 
on the collective motion. We consider the case where the coupling F is a one-body operator. 
We then need only the single-particle density matrix p of the intrinsic degrees of freedom to 
calculate the force. As we discuss below, temporal evolution of the single-particle density 
matrix exhibits a stochastic behavior. As a result, the intrinsic force has a fluctuating part 
on the top of its average value. Here, we find it more convenient to calculate the fluctuating 
part of force, which is determined by the fluctuating part of the single-particle density matrix 
Sp(t) = pit) — p[t), where the bar means taking ensemble average. Assuming fluctuations 
are small, 5 pit) is determined by a linearized transport equation [14] around the average 

Pit), 

thpp(t) - [h(t), 5p(t)\ - [5q(t)F, p(t)] = (3) 

where h(t) — h + q(t)F, h being the Hartree-Fock Hamiltonian of the separated nuclei and 
Sq(t) = q(t) — qit) denotes the fluctuation of the collective variable around its average value 
pit). The ensemble average value of the density matrix is determined by, 

ih^p(t)-[h(t),p(t)] = 0. (4) 

For simplicity of derivation, we neglect the collision term on the right hand side of eq.(3) and 
eq.(4) [15-17], however subsequently, we incorporate the damping width of single-particle 
states. Starting from an initial state 5p(s) at some time s, the formal solution of eq.(3) can 
be given as, 
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5p(t) = - % - I dt'5q(t') \G(t, t')FG ] (t, t'), pit)] + G(t, s)5p(s)G^(t, s) (5) 

it J s ^ ^ 

where the first term describes the effects of the perturbation during the time interval t— s with 
G(t, s) = exp —(i/h) J* dt'h(t') ~ exp —(i/h)(t — s)h(t) as the mean-field propagator, i.e. 
the propagator in the absence of thermal fluctuation, and the second term represent the 
propagation of the initial fluctuations Sp(s) of the intrinsic degrees of freedom during the 
time interval from s to t. Here, the initial time s does not represent the remote past, but 
rather it is sufficiently close to the time t, so that the time interval is much shorter than 
the relaxation time of the intrinsic degrees of freedom t — s <C r re j. Hence, the neglect of 
correlations due to collision term is justified in the description of eq.(5). In this case, the 
effect of correlations enters through the initial fluctuation term. Furthermore, we assume 
that the collective motion is sufficiently slow so that the intrinsic degrees of freedom is close 
to local equilibrium for each value of the collective variable. In order to evaluate the matrix 
elements of the fluctuating part of the density operator based on eq.(5), we approximate 
the average density matrix in terms of the instantaneous single-particle wave functions as 
P~(t) ~ \4 > i(t)) n i{4 > i(t)\i where we neglect the off diagonal elements. In this expression, 
the instantaneous wave functions (f>i(t) = (pi[q(t)] are determined from (h + qF)\(f>i(q)) = 
Ei(q)\4>i(q)) for each value of the collective variable q, and rij = 1/ [exp [(£j — ep)/T] + 1] 
denotes the Fermi-Dirac occupation factor at a temperature T. Employing the instantaneous 
representation, the matrix elements of fluctuations can be expressed as, 

SPij(t) = -^jyt'5q(t > )G ij (t : t')(i\F\j)(n j -n i ) + G ij (t : s)Sp ij (s) (6) 

with Gij(t, s) = exp[— i(t — s)(Ei — Sj)/h}. It is not possible to determine detailed structure 
of initial fluctuations of the intrinsic degrees of freedom. Therefore, it is plausible to assume 
that each matrix element of 5p(s) is a Gaussian random quantity with zero mean (i\Sp\j) = 
and a second moment. In accordance with the fluctuation-dissipation relation of the single- 
particle density matrix, we specify the second moment of Sp(s) as, 



{i\8p\j)(j'\6p\i') = 6n>8 jr - [m(l - rij) + nj(l - m)] . (7) 

In the special case of diagonal elements this formula gives the known result for fluctuations of 
occupation numbers, (i\5p\i)(i\8p\i) = n»(l — rij) [18]. Substituting eq.(5) into the right hand 
side of eq.(2), we find a generalized Langevin equation for the fluctuations of the collective 
variable, 

^5p(t) ± Mn 2 5q(t) = f dt'^t - t')5q(t') + £(t) (8) 
dt J s 

where the memory kernel in the retarded force and the random force term are given by 

7 ( f _f) = ^T,\(i\F\j)\ 2 Gji(t : t') Ml - nj ) -71,(1 -m)] (9) 

and 

ttt) = -Z(i\F\j)G ji (t,s)(j\5p\i). (10) 
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Using eq.(7), the auto-correlation function of the random force can be expressed as, 



mm = m\F\j)\ 2 Gij(t, mi - nj ) + n,a - ni )\ . (n) 

Dissipation and fluctuation aspects of dynamics are closely connected to each other, the 
similarity of expressions for the correlation function and the memory kernel reflects this fact. 
If the decay time of the memory kernel is sufficiently short, we can explicitly incorporate 
the memory effect into the retarded force in eq.(8). For evolution over a short time interval 
from i! to t, by neglecting the right hand side of eq.(8), we find the following relation, 

5q(t') « C(t - t')5q(t) - S(t - t')5p(t) . (12) 

For a parabolic potential well, propagators C(t — t') and S(t — t') are given by, 

C(t - t') = cos n(t - t') and S(t - t') = sin Q(t - t') . (13) 

On the other hand, for a parabolic potential barrier, these propagators are given by, 

C(t - t') = cosh n{t - t') and S(t - t') = sinh Q(t - t') . (14) 

The first term in eq.(12) involving 5q(t) introduces a shift in the curvature parameter of the 
potential. Here, we neglect this effect and substitute the second term on the right hand side 
of eq.(12) into the right hand side of eq.(8). Note that, since the fluctuations are linear, 
the same equation as eq.(8), but without the last term on the right hand side, holds for the 
average evolution by replacing 5p(t) and 5q(t) with p{t) and q(t), respectively. Therefore, 
combining the average evolution with the fluctuations, we obtain a generalized Langevin 
equation for the actual variables, p(t) = p(t) + 5p(t) and q(i) = q(i) + 5q(t), 

j t p(t) ± Mn 2 q(t) = -(3p(t) + t(t) . (15) 
Here, the reduced friction coefficient is given by, 

P= l f l ~ S ' dTE\(i\F\j)\ 2 e-i T ^- e ^n i {l-n i )S{r) + c.c.. (16) 
n Jo 

Substituting eqs.(13) and (14) for S(t), we find for the friction coefficient for a parabolic 
well, 



P = — Z\(i\F\j)\'' 



-%(t-8)(ej-£i-htl- 



Ej — Si — hQ — ir\ 

e -jr(t—s)(£j-£i+'hri-irf) _ 



Sj — £j + Ml — ir) 

and for a parabolic barrier, 



7ij(l - rij) + c.c. (17) 
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p 



2Mtt 



e _i( t _ S ) (£j ._ £i _j M -j ?7 ) _ ]_ 

Sj — Si — ihVL — if] 

g- jr{t—s)(ej—Ei+itel—vq) _ 

Sj — Si + — ii] 



rij(l — rij) + c.c. 



(18) 



In obtaining these results, we include the damping width rj of the single-particle states into 
the propagator in eq.(16) [19]. In further evaluation of the friction coefficients, we neglect 
the time-dependent terms in eqs.(17) and (18). The reason is the following: the dominant 
contributions to the friction coefficient arise from the coupling matrix element over an energy 
interval of the order of major shell spacing, Ej — = A m 10 MeV, which is much larger 
than typical values of collective frequency we consider here, HQ ~ 1.0 MeV. If the single- 
particle spectrum is sufficiently dense, the summations over the single particle states can 
be converted to energy integrals. As a result, exponential factors in eq.(17) and eq.(18) 
damp out over a time interval of the order of r = h/ A. Furthermore, in particular for low 
frequency collective motion, HQ < rj, these exponential factors damp out even over a shorter 
time scale as a result of damping of the single-particle states. Therefore, for a sufficiently 
long time interval, t — s ^> r , neglecting time dependent terms, we have for a parabolic 
well, 



(3 = mwm 



i 



V 



Mtt 

and for a parabolic barrier, 

p = mw\j)' n ' 



Ej — £j — Ml) + rf {Ej — Ei + Ml) + rf 



7lj(l - Uj) 



(19) 



{Ej - E t y + (m - v y {Ej - E,f + {mi + v y 



ni (l - rij) . (20) 



As seen from these results, for finite Q, the friction coefficient has different expressions 
around a well and a barrier. However, in the limit Q — > 0, it can be easily seen that 
these expressions become identical, known as the one-body friction formula. We call this 
limiting value as the classical friction coefficient and denote as @ . We introduce the Fourier 
transform of the correlation function of the random force [20], 



mm = j_ 



— -zu,(t-t ) cot h . 2D(u) 
x, 2tt 2T 2T v ; 



(21) 



where, 



D(u) = TE\(i\F\j)f 



CO 



rj 



(Ej - £i- hoo) + rf (Ej - Ei + huo) + rf 



nAl- rii 



(22) 



At low frequencies, D(oo) is just D(uo — > 0) = D = MT/3 the classical diffusion coefficient. 
On the other hand, the high frequency behavior is restricted by the magnitude of the cou- 
pling matrix elements. If the single particle spectrum is sufficiently dense, the magnitude of 
coupling matrix elements must decrease as a function of energy difference, mainly due to the 
mismatch of the overlap of the wave functions. We can represent this behavior by a Gaussian 
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or a Lorentzian function (i\F\j) 2 oc exp[— (ej — £i) 2 /2A 2 ] or oc 1/ [1 + (ej — £j) 2 /2A 2 ]. Fur- 
thermore, because of the Lorentzian factors in eq.(22), we can replace the energy difference 
Ej — Ei with the frequency Tiuu and approximately describe frequency dependence of diffusion 
coefficient as D{uS) = D exp [— (huj) 2 /2A 2 }, here we take the Gaussian for the frequency 
spectrum. As a result, the correlation function eq.(21) of the random force can be expressed 
as, 



(23) 



where 



x(t - f ) 



-°° duj 
=o 2tt' 



-tu(t-t) coth 



2T 



2T 



exp 



-(M 2 /2A J 



(24) 



The correlation function is characterized by two different parameters, the cut-off energy and 
temperature. Figure 1 shows the correlation function versus time at different temperatures 
T = 0.5 MeV, 1.0 MeV and 5.0 MeV. The results presented in this work are not very 
sensitive to the cut-off energy over a range of values A = 10 — 20 MeV. Therefore, in this 
figure and all other figures, we employ A = 15 MeV for the cut-off energy. At relatively high 
temperature hu <C 2T, (huj/2T) coth(huj /2T) ks 1, and the correlation function reduces to 
its classical form, 



, N f+°° duj 



e ^'-^exp 



2A 2 



1 



2ht ( 



exp 



-{t-t'Y/2r 2 



(25) 



For sufficiently short decay time To, it can be approximated by a delta function, Xo{t — t') — > 
5(t — t'), and as a result at high temperatures, the temporal evolution becomes Markovian 
and the random force acts like a white noise. However, as seen from the figure, in the 
quantal regime i.e. at low temperatures TiVL > 2T, we are faced with a stochastic evolution 
with a correlated noise [11,12]. 



III. ANALYSIS OF LANGEVIN EQUATION 

In order to obtain the joint probability distribution function P(q,p,t) of the collective 
variable and its conjugate momentum (q,p) by numerical simulation of the Langevin equa- 
tion, in general, we need to generate a sufficiently large ensemble of trajectories. Since, we 
have a correlated noise problem, we cannot use the standard methods [17,21,22] and we need 
to adopt suitable algorithms for numerical simulations [23]. However, in the situation that 
we consider here, the solution of the Langevin eq.(15) can be given analytically [13]. Since 
the equation is linear with a Gaussian random source, the probability distribution P(q,p,t) 
of collective variables is also Gaussian, which is determined by the mean values q(t), p(t) 
and the variances a qq (t) = 5q(t)5q(t), cr qp (t) = 5q(t)5p(t), <J pp {t) = Sp(t)Sp(t) of collective 
variables according to, 



2X 2 



((<? - (ifOqq + % -q)(P~ P)°q P + (P ~ Pf&pp) 



(26) 
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where X 2 = a qq a pp — a 2 p and is the inverse of the 2x2 matrix (cr^) with elements 
Cn — a qqi a \i — a qpi &i\ = °pg and °"22 = & PP - The mean values of collective variables q(t), 
p(t) are determined by the classical equations of motion, 

j t Q(t) = jjP(t) and j t p(t) ± MQ 2 q(t) = -fflt) . (27) 

Equations for variances are deduced from the Langevin equations for the fluctuating quan- 
tities Sq(t) = q(t) — q(t) and 5p(t) = p(t) — pit), 

j t Sq(t) = ±8p(t) and ^5p(t) ± Mn 2 5q(t) = -(35p(t) + t(t) . (28) 

Multiplying both sides of these equations by 5q(i), Spit) and performing ensemble averaging, 
we find 

j t <r qq (t) = ^a qp (t) (29) 
j t a qp (t) ± MQ 2 a qq (t) = ^a pp (t) - $a qp {t) + D qp (t) (30) 
jfrpit) ± 2Mn 2 a qp (t) = -2(3a pp (t) + 2D pp (t) (31) 



where D pp {t) = Sp(t)£(t) and D qp (t) = Sq(t)£(t) denotes the momentum and mixed diffusion 
coefficients, respectively. In order to evaluate diffusion coefficients, we need to calculate the 
dynamical fluctuations of collective variables in terms of the random force. This is carried 
out in Appendix A. Using the results for Sp(t) = J* * dt'Q(t — t')£(t') and from (A5) and (A6), 
diffusion coefficients can be expressed in terms of the correlation function of the random 
force as, 



D pp (t) = J o dt'Q(t - t') ■ £(t%(t) = 2D J o dsQ{s) ■ x (s) (32) 



and 



D qp (t) = f dt'S(t - t') ■ at')Z(t) = 2D f ds S(s) ■ x(s) . (33) 
Jo Jo 

In these expressions, the initial time is taken to be zero for convenience, t = 0, and the 
propagators Q(s) and S(s) associated with collective variables are given by (A7) and (A8) 
in Appendix A. At sufficiently high temperatures, correlation function x( s ) can be approx- 
imated by a delta function, and consequently, the momentum diffusion coefficient is time 
independent and takes its classical value, D pp = D and furthermore the mixed diffusion co- 
efficient vanishes, D qp = 0. The mixed diffusion coefficient is a genuine non-Markovian term, 
and it is absent in the Markovian limit. At low temperatures, due to the non-Markovian 
behavior of the correlation function, diffusion coefficients become time dependent and their 
magnitude are strongly modified by the quantum statistical fluctuations. We, also, note 
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that the modified frequency Vt = VL 2 + ((3 / 2) 2 enters in propagators Q(s) and S(s). For 
typical values of the frequency parameter HQ ~ 1.0 MeV of the conditional saddle and the 
magnitude of the reduced friction coefficient hf3/2 pa 1.7 MeV are comparable. As a result, 
the friction coefficient introduces a sizable modification of diffusion coefficients, which was 
not incorporated in the previous investigation [9,10]. Figures 2 and 3 show the diffusion 
coefficients in units of D , i.e., D pp / D and D qp /D , as a function of time for different values 
of temperature, T = 0.5 MeV, T = 1.0 MeV and T = 5.0 MeV. In order to illustrate the 
effect of friction on diffusion coefficients, we calculate diffusion coefficients by replacing Q 
with Q in the propagators Q(s) and S(s). Figures 4 and 5 compare two different values of 
diffusion coefficients D pp / D and D qp /D calculated with fl and Q as a function of time at 
temperature T = 1.0 MeV. The variances a qq , a qp and cr pp can be determined by solving 
the coupled differential equations (29), (30) and (31). However, it is much easier to obtain 
these variances directly from the Langevin equation (15) with the help of one-sided Fourier 
transform [18], as shown in Appendix A. 



IV. QUANTUM STATISTICAL EFFECTS ON DIFFUSION ALONG 
CONDITIONAL SADDLE TOWARDS FUSION 



In this section, we apply the generalized Langevin approach to investigate the influence 
of quantum-statistical fluctuations on diffusion along the fusion barrier, i.e. the formation 
probability Pf(t) of compound nucleus. When the conditional saddle, i.e. the inner fusion 
barrier, is approximately represented by an inverted parabola, the formation probability, 
i.e. the probability to cross the saddle point, can be calculated analytically in terms of 
distribution function of the elongation parameter q as [5,24,25], 

f°° 1 
Pf(t) = / dq 

Here, q(t) and cr qq {t) are the mean value and the variance of the elongation parameter, 
which are given by (A10) and (All) in Appendix A. In these expressions, (<? ?Po) are the 
mean values of elongation parameter and its conjugate momentum, and (a q0 , <r p0 ) are the 
associated variances at the initial configuration. As it is stated already, during the approach 
phase of the collision, system overcome the Coulomb barrier and some of the initial kinetic 
energy is dissipated into internal excitations and a sticking configuration is formed. The 
quantities (q ,p ) denote the average values of the elongation parameter and its momentum 
at the sticking configuration. In the second stage of the process, the shape of the system 
evolves from a sticking di-nuclear configuration towards formation of a spherical compound 
nucleus or re-separate again. The asymptotic value Pf(t — > oo) gives the transmission 
probability from a di-nuclear configuration to compound nucleus. In order to compare the 
results with our previous calculations [10], we consider collision of 48 C a and 238 U nuclei and 
adopt the same value for the reduced friction coefficient and the curvature parameter of the 
conditional saddle to be (3 = 5 x 10 21 s _1 and HQ = 1.0 MeV. We choose the initial position 
g to make the height of the conditional saddle to be 4.0 MeV and neglect the dispersion, 
i.e. (a q o,a p o), in the initial configuration. In the classical limit, the variance cr qq (t) of the 
elongation parameter has a analytical expression given by (A14), while in the quantum limit 



(g~g(*)) 
2a qq (t) 



-erfc 



2<7 M (t). 



(34) 
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it is given by (A18), and involves a one dimensional numerical integration over the frequency 
to. Figure 6 shows the formation probability Pf(t — > oo) of compound nucleus as a function 
of the initial kinetic energy K = p^/2M relative to the fusion barrier Vb at temperatures 
T = 0.5 MeV, T = 1.0 MeV and T = 5.0 MeV. These results, which are not very sensitive 
to the cut-off factor A, are presented for A = 15 MeV. Solid lines and dashed lines show 
the quantum and the classical calculations, respectively. At low temperatures, the quantum 
statistical fluctuations give rise to an enhancement of the formation probability, which is 
relevant to synthesis of superheavy elements by heavy-ion fusion reactions. The quantum 
enhancement is slightly less pronounced than that in the previous calculations [10]. The 
difference arises from the fact that in previous calculations the mixed diffusion coefficient 
D qp , which is a genuine non-Markovian term, is neglected and the momentum diffusion 
coefficient D pp is calculated with the unperturbed frequency Q, rather than Q. 



V. CONCLUSIONS 



Since several years, much effort has been made for synthesis of superheavy elements by 
heavy-ion fusion at near barrier energies. Even though, the reaction mechanism of heavy-ion 
fusion is not well understood, in most theoretical descriptions fusion is viewed as a diffu- 
sion process, which can be described by a Fokker-Planck approach or a stochastic Langevin 
approach [2-5,24,25]. In these descriptions the quantum statistical effects are ignored and 
a classical treatment is employed. In a recent work, we introduced a description based 
on a generalized Fokker-Planck approach, which incorporates quantum effects through non- 
Markovian transport coefficients [9,10]. In the present work, we follow a different description 
based on a generalized Langevin approach. The friction term and the random force both 
involve memory effects and they are related to each other in accordance with the fluctuation- 
dissipation theorem. As a result, the quantum statistical fluctuations are incorporated into 
the description. In principle, both approaches provide an equivalent description, however 
the Langevin approach has certain advantages in practical applications. In general, for a 
complex potential energy surface, it is much easier to carry out simulations of Langevin 
equation than solving the Fokker-Planck equation. Furthermore, diffusion coefficients pre- 
sented here are strongly modified by the friction mechanism, which was not considered in 
the Fokker-Planck approach [10]. In this work, we consider a simple model in which the 
fusion barrier is represented by an inverted parabola. In this case, the joint distribution 
function of the collective variable and the conjugate momentum becomes a Gaussian, and 
the fusion probability can be given in an analytical form. Calculations illustrate that the 
quantum statistical effects enhance the fusion probability at low temperatures. However, 
the enhancement is somewhat less pronounced than that reported in the previous investi- 
gation [10]. For a realistic potential energy surface, analytical solution is not possible and 
distribution function of collective variables should be constructed by generating sufficient 
number of events of the generalized Langevin equation. Since, the random force is not a 
white noise, we cannot use standard methods for numerical simulations [21,22]. Therefore, 
it is necessary to develop suitable algorithms for simulations of the Langevin equation with 
a correlated random force [23] . 
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APPENDIX 

In this Appendix, we analyze solutions of the Langevin eq.(15) together with dq/dt = 
p/M by employing one-sided Fourier transform [18]. After performing the Fourier transform, 
we obtain 

—go — iujq{uj) — an d — Po — iujpiuj) ± MVL 2 q{uj) = —j3p(u) + £(cu) (Al) 

where (qo,Po) are the initial conditions, qiu) = / °° dtexp(iut)q(t) is the one-sided Fourier 
transform of the coordinate. The p(oo) and C,(uj) are similarly defined. Combining two 
equations in (Al), the Fourier transforms of the collective variables are given by, 

/ \ uj + i(3 1 p + £(uj) 

9M=29 V T ^ + ^ "M ^ T ^ + ^ (A2) 

p(w) = ^wT^p + ' (A3) 

Time dependence of the collective variables are found by the inverse Fourier transformation, 

/+°° diiJ r+oo (J^fjj 

— exp(—iut)q(u) and p{t) = / — exp(— iut)p{u) . (A4) 

Integration over uj in these expressions can be carried out with the help of Cauchy theorem. 
Here, we give the results for a parabolic potential barrier, 

q{t) = q R(t) + p S(t) + f dt'Sit - 0£(0 (A5) 

Jo 

and 

p(t) = q (Mnf S(t) + p Q(t) + f dt'Q(t - 0£(0 (A6) 

Jo 

where Q(t), S(t) and R(t) are given by, 

Q(t) = exp ^~^tj ^coshHt - ^= sinhfi^ (A7) 

s{t) = m exp {~^) sinhUt (A8) 
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and 



R(t) = exp ( — ( cosh ilt + ^Lsinhfit 



(A9) 



where Q — y f2 2 + (f3/2) 2 . The solutions can be given in a similar manner for a parabolic 
potential well. The mean values of collective variables are obtained by taking the ensemble 
average of (A5) and (A6), 



q(t) = q R(t) + p S(t) and p(t) = q (Mfi) 2 S(t) + p Q(t) . 
The variances are given by, 

<r qq (t) = (r q0 R 2 (t) + a p0 S 2 (t) + f ds f ds'S(s)S(s')2D oX (s - s') 

JO JO 



(A10) 



(All) 



(Tgpit) = a q0 (MQ) 2 R(t)S(t) + a p0 S(t)Q(t) + jf* rfs jf* ds / S , (s)Q(s')2A,x(s - s') (A12) 
and 

MO = (Mfi) 4 S 2 (t) + a p0 g 2 (t) + ds J* ds'Q(s)Q(s')2D oX (s - s') . (A13) 

In these expressions, first two terms describe propagation of the initial fluctuations of the 
coordinate and momentum distributions a q o, cr p o and the last term arises from dynamical 
fluctuations generated by the random force. 

For calculation of the formation probability of compound nucleus, we need an explicit 
expression for variance cr gq (t) of the collective variable. In the Markovian limit, using the 
fact that the correlation function behaves like a delta function, X (s ~ s ') ~^ <K S ~~ s ')> we 
obtain the known analytical result for the dynamical part of the cr qq (t) [5], 

cr* (£) = jT ds J* ds'S(s)S(s')2D oX (s - s') -> J*dsJ*ds'S(s)S(s')2D 5(s - s') 



T 



-^(sinhfit) 2 + -^=(sinh2fit) - e +/3t + 1 
2Vt 2\l 



(A14) 



In the classical limit, the dynamical part of the cr qp (t) and (J pp (t) are similarly given by, 



= f ds J*ds'S(s)Q(s')2D oX (s - s') -> jf* ds J* ds'S(s)Q(s')2D 5(s - s') 



^e'^smhUt) 2 



(A15) 



and 



= ds ds'Q{s)Q{s')2D QX {s - s') -> J* ds J* ds , Q(s)Q(s')2D 5(s - s') 



MTte-P* sinhfit 



cosh fit sinh fit 

2fi 



+ MT (l — e~^) . 



(A16) 
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For quantal calculations, introducing the Fourier transform of the correlation function, 

f+°° duo 



(A17) 



we can express the dynamical part of the variance in terms of a one-dimensional numerical 
integration over the frequency u as, 



-oo A, \ _ o 

— S t (u>) x(oo)2D 

x> 27T 



(A18) 



where S t (uj) = J * dsS(s)e luJs . Dynamical parts of variances <J qp (t) and <y pp (t) can be evalu- 
ated in a similar manner to give, 



^ — 5 t (o;)g t *(a;)xH2L> and = jf^ — xM2D (A19) 

where Q t (u) = Jq dsQ(s) e - iu)S . 
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FIGURES 
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1. The correlation function is plotted versus time for A = 15 MeV. 
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FIG. 3. The mixed diffusion coefficient, in units of classical diffusion coefficient, is 
plotted versus time for A = 15 MeV. 
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FIG. 4. The momentum diffusion coefficient, in units of classical diffusion coefficient, is plotted 
versus time for A = 15 MeV and T=1.0 MeV. The cases with and without friction are compared. 
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FIG. 5. The mixed diffusion coefficient, in units of classical diffusion coefficient, is plotted 
versus time for A = 15 MeV and T=1.0 MeV. The cases with and without friction are compared. 




FIG. 6. The formation probability is plotted versus initial kinetic energy minus barrier height. 
The result of classical diffusion approach is compared with that of quantum diffusion approach for 
A = 15 MeV. 



16 



